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We propose a scheme to demonstrate fractional statistics of anyons in an exactly solvable lattice 
model proposed by Kitaev that involves four-body interactions. The required many-body ground 
state, as well as the anyon excitations and their braiding operations, can be conveniently realized 
through dynamic laser manipulation of cold atoms in an optical lattice. Due to the perfect local- 
ization of anyons in this model, we show that a quantum circuit with only six qubits is enough for 
demonstration of the basic braiding statistics of anyons. This opens up the immediate possibility of 
proof-of-principle experiments with trapped ions, photons, or nuclear magnetic resonance systems. 



PACS numbers: 

Anyons, as exotic quasiparticles living in two dimen- 
sions with fractional statistics [![, have attracted strong 
interest over the past two decades. The excitations in 
fractional quantum Hall systems have been predicted to 
be anyons [2j , and fractional charges have been confirmed 
with some experimental evidence Q. However, a direct 
observation of fractional statistics associated with anyon 
braiding is hard in this system and has recently attracted 
intriguing theoretical proposals 0]. Other systems with 
anyon excitations have also been proposed. In partic- 
ular, in the context of topological quantum computa- 
tion, Kitaev described two exactly solvable theoretical 
models 0,0] which support anyons. An implementation 
scheme has been proposed for the second model Hamilto- 
nian 0, Hj], which involves only two-body interactions. A 
detection method of anyons associated with this imple- 
mentation scheme has been proposed recently Q. This 
implementation scheme with an optical lattice requires 
achievement of very low temperature 0] , which unfortu- 
nately is still somewhat beyond the current experimental 
capabilities. 

In this paper, we propose a method to realize anyons 
and to observe their braiding statistics associated with 
the first Kitaev model that is based on four-body inter- 
actions. Four-body interactions are notoriously difficult 
to generate experimentally in a controllable fashion. Our 
key idea here is to generate dynamically the ground state 
and the excitations of this model Hamiltonian instead of 
direct ground-state cooling. Note that the anyon proper- 
ties are associated with the underlying many-body entan- 
gled states. If we generate exactly the same ground and 
the excited states, we should be able to observe the same 
property. We propose several systems for implementa- 
tion of this idea. First, we show that laser manipulation 
of ultra-cold atoms in an optical lattice provides a nat- 
ural system to realize a large scale ground state of this 
model Hamiltonian. The braiding statics of anyons can 
then be observed with a series of single-bit operations. 
Then, we show that the systems of trapped ions [icj |. 
photons [Til ], or nuclear magnetic resonance 12 1, pro- 



vide a ready platform to realize a small scale system for 
proof-of-principle demonstration of the anyon braiding 
statistics. The anyons are perfectly localized quasiparti- 
cles in this model Hamiltonian, which means we do not 
need a large system for implementing their braiding op- 
erations. We show that as little as six qubits are enough 
for demonstration of the basic braiding statistics; and 
with nine qubits, one can further demonstrate robust- 
ness of the braiding operation to certain variations of the 
braiding path. 

The first Kitaev model is a spin Hamiltonian for a two- 
dimensional square lattice [H[ (more generally, we can 
extend the model from a square lattice to any planar 
graph). One associates each edge of the lattice with a 
qubit (a spin-1/2 particle). The model Hamiltonian is 
given by 
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where A v 
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is defined for each vertex v as A v — 
and Bf is defined for each face / as 

/ = Uj=boundary(f) Z J ( See Fi §- 1A for an illustration). 
The operators X and Z denote the standard Paul ma- 
trix o x and <7z, respectively, and A v and Bf are called the 
stabilizer operators in the context of quantum error cor- 
rection. In a square lattice, each term of the Hamiltonian 
H represents four-body interaction of local qubits. One 
can easily check that all the terms of this Hamiltonian 
commute, so it is straightforward to get the ground state 
\ip) of the Hamiltonian H , which is given by A v \<p) = \<p) 
and Bf\tp) = \(p) for all vertexes and faces. The state \<p) 



is highly entangled [13 1. 

A quasiparticle is generated on the vertex vq (or face 
fo) if A Vo (or Bf ), acting on the excited state \<p e ), 
yields an eigenvalue —1 instead of +1 for the ground 
state. The quasi-particles in these two cases are called e- 
particles (vertex) or m-particles (face), respectively. One 
can check that the e-particles and m-particles by them- 
selves are bosons, but the mutual statistics between the e 
and m particles become the one for 1/2-anyons 0, as we 
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get a phase flip e l7r if we move the e(m)-particle around 
the m(e)-particle along an arbitrary loop. 

It is hard to directly generate the interactions repre- 
sented by the Hamiltonian H and to cool the system to 
its ground state \ip). However, we note that the braid- 
ing statistics of anyons are directly associated with the 
entanglement properties of the underlying ground and 
excited states, whereas the Hamiltonian only plays an 
implicit role for nailing down the corresponding state. 
As long as we can create the state \ip) and generate the 
e and m-particles above this state, we should be able 
to demonstrate the fractional braiding statistics between 
these quasiparticles. Our task reduces to how to cre- 
ate the state \tp) and how to demonstrate the fractional 
statistics of quasiparticles excited from the state \ip). In 
the following, we propose two approaches, targeted at 
a large scale implementation and a small scale proof-of- 
principle demonstration, respectively. 

Physically, the simplest method to generate a large 
scale state \<p) is to start from a two-dimensional (2D) 
cluster state. The cluster and the graph states, intro- 
duced in [3, EH , are defined as the co-eigenstate of a set 
of commuting stabilizer operators Si (with +1 eigenval- 
ues). The 2D cluster state is associated with qubits on 
vertexes of a 2D square lattice, where a stabilizer opera- 
tor is introduced for each vertex as Si — X; ®jejv(,*) 
Therein, N(i) denotes the set of nearest neighbors of the 
vertex i. A 2D cluster state can be conveniently gen- 
erated with ultra-cold atoms in a square optical lattice. 
Each atom is effectively two-level, which defines a qubit 
(or a spin- 1/2 particle). One starts with all the atoms in 
equal superpositions of |0) and |1) states (co-eigenstates 
of Xi). Then, with control of the optical lattice potential, 
one turns on spin-dependent collisions [lj| or tunneling 
for a fixed amount of time (both give rise to an ef- 
fective Ising interaction) to get the 2D cluster state 14 1 . 
Starting from a 2D cluster state, single-bit measurements 
of half of the qubits in the bases X and Z respectively 
will yield the desired state \tp) [l7j |. The measurement 
pattern is illustrated in Fig. IB. One can check these 
measurements transfer the set of the stabilizer operators 
from {Si} to {±A U , ± Bf} for the remaining qubits, as 
it is required for the ground state of the Hamiltonian 
(1). The sign factors ±1 in the stabilizer generators de- 
pend on the measurement outcomes from the measured 
qubits. All these sign factors can be readjusted to +1 by 
subsequent local Pauli operations. 

With creation of the state \ip), one can demonstrate 
the fractional statistical phase of the 1/2-anyons through 
a Ramsey type interference experiment [{J fig ]- First, 
by applying a single-bit rotation Z to one of the edge 
qubits, one creates a pair of e-particles on the neighbor- 
ing vertexes, as represented by an excited state \ip e ) of 
the Hamiltonian (1). Then, with half of the rotation 
y/~Z, we get the superposition Qtp) + \ip e }) /V%- Simi- 
larly, with a single-bit rotation X on another edge qubit, 




FIG. 1: (Color Online) (A)An illustration of the Kitaev model 
Hamiltonian; (B) The measurement pattern to prepare the 
many-body ground state \<p) from a two-dimensional cluster 
state, where X and Z denote respectivel y a single-bit mea- 
surement in the X or Z basis, (see also [13]); (C) An illus- 
tration of braiding of anyons, where a pair of e-particles (m- 
particles) are connected respectively by the ci (or oi) lines, 
and c denotes the braiding loop. The X or Z besides the edge 
qubits denote the corresponding single-bit operations for cre- 
ation and braiding of the anyons. See the text for more de- 
tailed description. 



we create a pair of m-particles on the neighboring faces. 
Then, one can move one of the m-particles around one 
of the e-particles along a loop through successive X ro- 
tations acting on the qubits in the loop (see Fig. 1C). 
After fusion of the m-particles (through another X ro- 
tation), the underlying state becomes (\(p) — \<p e )) /V% 
due to the fractional phase 7r acquired from braiding of 
the e and m-particles. This phase flip can be unambigu- 
ously detected. For instance, with another \f~Z operation, 
(\<p) — \<p e }) /v2 goes to \<p) (whereas without the statis- 
tics phase 7T, + \<p e )) /v2 would go to \<p e ))- The 
states \<p) and \ip e ) can be distinguished by measuring 
the relevant stabilizer operators {A v } (measurement of 
the stabilizer operators only requires single-bit detection 
together with classical correlation of the measurement 
outcomes). 

In the above implementation scheme, except for the 
initial step of preparation of the cluster state with spin- 
dependent atomic collisions, all the following steps are 
based on single-bit unitary operations or measurements. 
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With cold atoms in an optical lattice, one can indeed re- 
alize a large-scale system with millions of qubits (l9j . 
and thus can well demonstrate the robustness of the 
braiding operation: the statistics phase it depends only 
on the topology of the loop and not on the detailed 
path. The only remaining experimental challenge in this 
scheme is the requirement of addressing of individual 
qubits for single-bit operations. The experiments on ul- 
tracold atoms are making progress towards the ability of 
single-bit addressing. On the other hand, there are also 
other experimental systems which start from the bottom 
up and have demonstrated the ability to fully control a 
small number of qubits. Along with that line, in the fol- 
lowing we propose a different implementation scheme to 
provide proof-of-principle demonstration of anyons and 
their braiding statistics in small systems. This scheme is 
well within reach of the current experimental technology. 

It is interesting to ask what the smallest system size 
is for demonstration of the braiding statistics of anyons. 
The system size has to be significantly larger than the 
size of anyons for the braiding operation. For any systems 
with anyons, the quasiparticlcs are necessarily localized 
in space (otherwise it is impossible to define the braid- 
ing). For the first Kitaev model, the quasiparticles are 
particularly well localized: the e-particle is on a single 
vertex and the m-particle is on a single face. Due to this 
perfect localization, we do not need a large system for 
a basic braiding operation. It turns out that the graph 
shown in Fig. 2A is the smallest system for implementa- 
tion of the anyon braiding operation, which just requires 
manipulation of six qubits. 

The graph in Fig. 2A corresponds to the ground state 
of the Hamiltonian H§ = —A\—A 2 — B\—B 2 — Bj ) — B4oi 
six edge qubits, where A\ = XiX 2 X 3 , A 2 — A3A4A5A6, 
B\ = Z1Z3Z4, B 2 = Z 2 Z 3 Z 5 , B 3 — Z^Zq, B4 — Z 5 Z 6 . 
One can check that its ground state \tp) e is equivalent 
under local single-bit operations to a graph state shown 
in Fig. 2B (see Ref. [l§] for definition of the graph state). 
One can use five controlled phase flip (CPF) gates to pre- 
pare the graph state of Fig. 2B and thus also \ip) 6 . Figure 
2C shows the detailed preparation circuit for the state 
|<^) 6 which involves five CPF gates as well as a few Hard- 
mard (H) gates. Then, as also shown in Fig. 2C, a v^3 
operation on the qubit 3 generates the superposition state 
(|<^) 6 + \(p e )a) /V2, where |(/J e ) 6 has a pair of e-particles 
on two vertexes of the edge 3. Another X± operation on 
the qubit 4 further generates a pair of m-particles on two 
neighboring faces of the edge 4. With four X operations 
on the qubits 6-5-3-4, we achieve the braiding by moving 
an m-particle around a e-particle along the loop shown 
in Fig. 2A. Finally, we fuse (annihilate) all the quasi- 
particles (with another \/~Z$ and X4 operation), and the 
resultant state is detected with single-bit measurements 
(in X or Z basis) of the six stabilizer operators. 

In the above six-qubit minimum implementation, the 
braiding loop is unique and there is no way to show the 
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FIG. 2: (Color Online) (a) A planar graph representation of 
the Hamiltonian with a minimum number of qubits that sup- 
port braiding of anyons. (B) The graph state that is equiva- 
lent under local unitary operations to the ground state of the 
Hamiltonian in Fig. 2A. (C) The detailed quantum circuit for 
generating the ground state \ f} 6 , for creation, braiding, and 
annihilation of the anyons, and for the final state detection. 



robustness of the topological braiding operation by mov- 
ing the anyons along different paths. To have demon- 
stration of some robustness of the anyon operation to 
different braiding paths, we give another implementation 
in Fig. 3 which uses nine qubits. The ground state of this 
nine-qubit Hamiltonian represented by Fig. 3A is locally 
equivalent to graph state shown in Fig. 3B, and the de- 
tailed implementation circuit is given in Fig. 3C. In this 
plane geometry (see Fig. 3A), there are several different 
loops for the quasiparticles. The anyonic operation only 
depends on the topological character of the loop, i.e., 
whether the different types of quasiparticle (e and m) 
"braid" with each other. For instance, both the loops 6- 
5-3-4 and 9-8-5-3-4-7 of the m-particles (the dashed lines 
in Fig. 3A) give the same result of the statistical phase, 
while the loop 9-8-6-7 (the solid line in Fig 3A) has no ef- 
fect since it does not enclose a e-particle. Braiding along 
different loops can be implemented with single-bit X op- 
erations on different sets of individual qubits. See Fig. 3C 
for an example circuit implementing braiding of a m- 
particle with an e-particle. The m-particle is taken along 
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FIG. 3: (Color Online) (A) A planar graph representation of 
the Hamiltonian with 9 qubits that support different types of 
loops for anyons. The dashed line loop denotes an effective 
braiding while the solid line loop is not. (B) The graph state 
locally equivalent to the ground state of the Hamiltonian in 
Fi.g 3A. (C)The detailed quantum circuit for realization of the 
ground state and the associated anyon braiding (the braiding 
corresponds to the dashed line loop in Fig. 3A). 



the loop 9-8-5-3-4-7. The enclosed e-particle resides on 
the vertex between edges 3,4,5,6. 

In the context of quantum computing, entangling gates 
of six to eight qubits have been demonstrated in with 
trapped ions 1 1 Oil . six-photon entanglement in the linear 
optics system and manipulation of a dozen of nu- 
clear spins in the nuclear magnetic resonance system [l2| . 
With such experimental capabilities, the above six-bit 
and nine-bit implementation schemes can be expected to 
be realized in the very near future. This kind of proof-of- 
principle demonstration of anyons in small and relatively 
simple systems will represent an important step towards 
the long pursued goal to demonstrate fractional statis- 
tics of quasiparticles in a macroscopic material. Such 
abilities, properly extended to large systems, will also be 



critical for future implementation of fault-tolerant topo- 
logical quantum computation. 

In summary, we have proposed a method to demon- 
strate fractional braiding statistics of anyons in an ex- 
actly solvable spin model. Two types of implementation 
schemes are described, respectively targeted at a large- 
scale realization with cold atoms in an optical lattice and 
a small-scale realization with some available qubit sys- 
tems, such as trapped ions, photons, or nuclear spins in 
liquids of molecules. The proposed schemes open the 
prospect of experimental demonstration of anyons in the 
near future. 
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